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Abstract 

We present a study of the results obtained by combining LO partonic matrix elements 
with either LO or NLO partons distributions. These are compared to the best prediction 
using NLO for both matrix elements and parton distributions. The aim is to determine 
which parton distributions are most appropriate to use in those cases where only LO 
matrix elements are available, e.g. as in many Monte Carlo generators. Both LO and 
NLO parton distributions have flaws, sometimes serious, for some processes, so a modihed 
optimal LO set is suggested. We investigate a wide variety of process, and the new modified 
LO* pdf works at least as well as, and often better than, both LO and NLO pdfs in nearly 
all cases. The LO* pdf set is now available in the LHAPDF package [1]. 



1 Introduction 



It has long been known that for certain regions of x there can be big differences between 
parton distributions extracted at different orders of perturbative QCD (see for example [2] for 
a discussion). Moreover these differences between LO, NLO, and NNLO partons are much more 
significant than those between parton distributions extracted by various groups. This happens 
due to important missing higher order corrections in the splitting functions which determine 
the parton evolution, and in the cross-sections which govern their extraction by comparison to 
experimental data (mainly from structure functions in deep inelastic scattering). In particular, 
use of parton distributions of the wrong order can lead to incorrect conclusions on the size of 

^ Royal Society University Research Fellow 
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Figure 1: A comparison of various gluon distributions at NLO and a gluon distribution at LO. 
There is significant variation in the NLO distributions, but the LO distribution is a quahtatively 
different shape from all of them. 
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the small-x gluon, shown in Fig. 1, and consequently on the importance of corrections due to 
shadowing, saturation etc. at small x. 

This large change in the size of parton distributions as one changes perturbative order means 
that the combination of the respective order of the parton distribution function (pdf ) and the 
accompanying matrix element is an important issue when calculating the cross-section for any 
physical process [3]. Traditionally, LO pdfs are thought to be the best choice for use with LO 
matrix elements, such as those available in many LO Monte Carlo programs, though it has been 
recognised that all such results should be treated with care. However, recently an alternative 
viewpoint has appeared, and it has been suggested that NLO pdfs may be more appropriate [4]. 
The argument is that NLO cross-section corrections are often small, and the main change in the 
total cross-section in going from LO to NLO is due to the pdfs. Indeed, there has already been 
an investigation of the use of NLO pdfs for the underlying event [5] at the Tevatron. There is a 
big difference in the results when using CTEQ6L and CTEQ6.1M pdfs [6] due to the changes 
in the gluon, though agreement can be reached by significant retuning. However, this retuning 
will potentially affect predictions for other quantities. 

In this article we investigate the differences in predictions obtained for a variety of physical 
quantities combining different pdfs with LO matrix elements. In each case we assume NLO pdfs 
combined with NLO matrix elements represent the best prediction - the truth? The different 
processes chosen rely on a variety of input parton distributions and on various ranges of x and 

in order for us to try to make our conclusions as all-encompassing as possible. We interpret 
the features of the results noting that there are significant faults if one uses exclusively either 
LO or NLO pdfs. We hence attempt to minimise this problem, and investigate how a best set 
of pdfs for use with LO matrix elements may be obtained. 

2 Parton Distributions at Different Orders 

Let us briefly explain the reasons for the origins of the differences between the parton distribu- 
tions at different perturbative order. The qualitative difference between the LO and NLO gluon 
distribution is shown in Fig. 1. Clearly the LO gluon is much larger at small x than any NLO 
gluon at = 5GeV^. The evolution of the gluon at LO and NLO is quite similar, so at larger 
the relative difference is smaller, but always remains significant. At = 1 — 2GeV^ the 
difference is even more marked, with NLO gluon distributions becoming valence-like, or even 
negative at very small x, while the LO distribution remains large. This difference in the gluon 
distributions is a consequence of quark evolution, rather than gluon evolution. The small-x 
gluon is determined by dF2/ dlnQ'^ , which is directly related to the evolution of the quark 
distributions. The quark-gluon splitting function Pqg is finite at small x at LO, but develops a 
small-x divergence at NLO (and further ln(l/a;) enhancements at higher orders), so the small 
X gluon needs to be much bigger at LO in order to fit structure function evolution. 

^Of course, the full NNLO calculation is known for some processes, but we will avoid this further complication 
in this article. 
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Figure 2: The comparison of the up quark parton distribution at LO with that at NLO in 
the MS factorization scheme (left) and the comparison between the NLO MS-scheme and DIS- 
scheme up quark distribution (right), both at at = 6400GeV^. The upper plots show the 
uncertainty bands of the NLO MS distribution and the lower plots the percentage uncertainty. 
In both cases this is compared to the central LO or DIS scheme distribution. 
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There are also significant differences between tfie LO and NLO quark distributions, as sfiown 
in tfie left of Fig. 2. Before analysing these in detail it is important to realise that beyond LO the 
parton distributions are factorization scheme dependent. We conventionally take NLO to mean 
"NLO in the MS factorization scheme" , and indeed this is the NLO distribution used in the left 
of Fig. 2. In order to illustrate the difference between orders we also show the comparison of the 
MS-scheme distribution at NLO and the DIS-scheme distribution in the right of Fig. 2. In the 
DIS-scheme the NLO coefficient functions relating parton distributions to structure functions 
are zero. Hence, if the LO parton distributions and NLO distributions produced identical 
structure functions, the right and left plots in Fig. 2 should be the same. Indeed, there are 
obvious similarities, both LO and NLO DIS-scheme quark distributions being larger than NLO 
MS-scheme distributions at very high x and very low x. This is due to the effect of the NLO 
coefficient functions in the MS-scheme. Most particularly the quark coefficient functions for 
structure functions in the MS scheme have ln(l — x) enhancements at higher perturbative order, 
and the high-x quarks are smaller as the order increases. However, the LO quark distribution 
is depleted compared to the NLO DIS-scheme quark distribution for 0.001 < a; < 0.1. And 
indeed, the LO fit to structure function data is poor in this region. 

Hence, while LO parton distributions are more similar to NLO DIS-scheme distributions 
than NLO MS-scheme distributions (the gluons in DIS scheme and MS scheme are similar 
except at very high x), there are some qualitative differences between LO and NLO in general. 
The LO gluon is much bigger at small x, and compared to MS scheme the LO valence quarks 
are much bigger at high x. This is then accompanied by a significant depletion of the quark 
distribution for x in the region of 0.01, despite the fact that this leads to a poor fit to data. 

Let us examine the consequence of these differences when calculating physical cross-sections. 
In Fig. 3 we illustrate the ratio of the cross-sections produced for Drell-Yan production of a 
boson with invariant mass SOGeV at Tevatron energies using NLO MS-scheme parton distri- 
butions and NLO MS-scheme matrix elements (taken as the reference) [7], NLO MS-scheme 
parton distributions and LO matrix elements and LO parton distributions and LO matrix 
elements.^ In this case for central rapidity, |/ = 0, both partons have x = 0.04, and the larger 
and smaller x, Xi and X2 are 0.04exp{y) and 0.04exp(— y) respectively. In this case it is the 
quark distributions which are probed. 

Clearly we are generally nearer to the truth with the LO matrix element and NLO parton 
distribution [8] than with the LO matrix element and LO parton distributions [9]. However, 
this is always too small, since the NLO correction to the matrix element is positive. It is 
relatively y-independent, but there is the beginning of extra suppression at the highest y for the 
PDF[NLO]-ME[LO] calculation due to an additional enhancement in the NLO matrix element 
at high y, similar to the ln(l — x) enhancement for structure function coefficient functions. 
The depletion of the LO quark distributions for x ~ 0.04 leads to the extra suppression in the 
PDF[LO]-ME[LO] calculation. However, when probing the high-x quarks, the increase in the 

^Sincc NLO matrix elements are most readily available in MS seheme, we will take this as the default, and 
henceforth NLO is intended to mean NLO in MS scheme unless explicitly stated otherwise. 
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Drell-Yan Cross-section at Tevatron for 80 GeV with Different Orders 
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Figure 3: Comparison of boson production at the Tevatron using combinations of different 
orders of parton distributions and fiard partonic cross-sections. 
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LO parton compensates for the increase in the NLO matrix element, and for y > 2 this gives 
the more accurate rcsTilt. However, overall the shape as a function of y is much worse using 
the LO parton distributions than the NLO distributions. 

We consider the same comparison, but for boson production at the LHC rather than the 
Tevatron, in Fig. 4. In this case central rapidity corresponds to a; = 0.006, and we have to 
go to higher rapidity than at the Tevatron to be sensitive to large- ,t effects, but have much 
more sensitivity to small x. Again the PDF[NLO]-ME[LO] calculation is suppressed compared 
to the truth, but this time there is less change in shape at high rapidity since even at y = 3 
{xi = 0.11, a;2 = 0.0003) we are barely reaching the regime of enhanced NLO partonic cross- 
sections. As for the Tevatron, the PDF[LO]-ME[LO] result is even more suppressed, due to the 
depletion of quarks, though this begins to reduce at the highest y, mainly due to increases in 
small-x quarks at LO. The general conclusion is the same as for the Tevatron - the NLO pdfs 
provide a better normalization and a better shape. 

This single quantity suggests that the opinion in [4] is correct, i.e. that NLO pdfs are more 
appropriate if one only has LO matrix elements. However, this is a very processes dependent 
statement. In order to demonstrate this we consider a quantity sensitive to the small- ,t ghion 
distribution, where the difference between the LO and NLO gluon distribution is due to missing 
higher order perturbative corrections. Let us consider the production of charm in DIS, i.e. 
F2^{x, Q^) where we take all charm as produced in the final state, i.e. we work in the fixed- 
ffavour number scheme (FFNS). In this case the NLO coefficient function, C^Jg'^\x,Q'^ ,m1) 
contains a divergence at small x not present at LO, in the same way that the quark-gluon 
sphtting function does, the latter being responsible for the large difference between the LO 
and NLO gluon distributions at small x. In the left of Fig. 5 we see the large effect of the 
NLO coefficient functions [10] . When using NLO partons the LO matrix element result is well 
below the truth at low scales, and the shape is totally wrong. The structure function follows 
the shape of the gluon distribution at NLO, whereas the small-a; divergence in the coefficient 
function provides a charm structure function which rises at small x in the same way that the 
NLO correction to Pqg provides a total dF2 / d\n{Q'^) which rises at small x. In the right of 
Fig. 5 we see the result using the LO pdfs. In this case F2''{x, Q"^) rises at small x even with 
the LO coefficient function, due to the rising gluon. Using the NLO coefficient function and 
the LO pdfs effectively double counts the small-x divergence, and the result is much too large 
and steep. All four results are shown together in Fig. 6. While the LO pdfs combined with LO 
coefficient functions is not a perfect match to the truth - after all the small- x divergences are 
not exactly the same in matrix element and splitting function - it is clearly far better than the 
other two approaches which either double count the small-a; divergences, or include neither. In 
particular, in this case the NLO pdfs together with the LO matrix elements fail badly, providing 
a clear contradiction to the view in [4]. 

Hence, from these two simple examples alone we can make the following conclusions. Some- 
times it is better to use NLO pdfs if only LO matrix elements are known, and one can encounter 
significant problems with both normalization and shape if LO pdfs are used. However, one can 
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Drell-Yan Cross-section at LHC for 80 GeV with Different Orders 
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Figure 4: The same as Fig. 3, but at LHC energies. 
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Figure 5: A comparison of charm production at HERA using LO and NLO coefficient functions 
and NLO parton distributions (left), and using LO parton distribution functions (right). 
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be completely wrong, particularly at small by using NLO partons with LO matrix elements, 
due to effective zero-counting of small-x divergences. One could finish here and conclude that 
it is necessary to consider on a case by case basis, and in any analysis which uses LO Monte 
Carlo generators sometimes use LO pdfs and sometimes NLO pdfs, depending on the particular 
process under consideration. However, we will continue our investigation rather than accepting 
this pessimistic conclusion and attempt to discover whether there are some optimal partons 
which have the desirable features of both LO and NLO pdfs. 

3 Improved LO Parton Distributions 

In order to make progress in finding some optimal set of pdfs for use with LO matrix elements we 
need to understand fully the differences between LO and NLO partons. As seen by comparison 
with the NLO DIS-scheme partons in Fig. 2, part of the dip in LO quarks compared to NLO 
is due to extra coefficient function contribution at NLO (particularly for x ~ 0.1), but it is 
mostly a problem at LO - the size and extent of the depletion from x — 0.1 — 0.001 is reflected 
in the fit quahty. At LO, compared to NLO (and higher orders), missing terms in ln(l — x) 
and ln(l/a;) in coefficient functions and/or evolution lead to a gluon which is much bigger as 
a; — > and valence quarks which are much larger as a; 1 in order to compensate. From the 
momentum sum rule there are then not enough partons to go around, hence the depletion in 
the quark distributions at moderate to small x. 

This depletion leads to a bad global fit at LO, particularly for HERA structure function 
data which are very sensitive to quark distributions at a; ~ 0.01. In practice the lack of partons 
at LO is partially compensated by a LO extraction of as{M'^) ~ 0.130, i.e. the very large 
coupling constant leads to quicker evolution, making up for the higher order corrections and/or 
smaller parton distributions in some regions. This results in one obvious modification. It is 
helpful to use the NLO definition of the coupling constant in a LO fit to parton distributions. 
Because of quicker running at NLO, LO and NLO couplings with the same value of 05 (M|) 
become very different at lower scales where DIS data exist, as shown in Fig. 7. Alternatively, 
near = IGeV^ the NLO coupling with as{M'^) = 0.120 is similar to the LO coupling with 
as{M'^) = 0.130. Consequently, the NLO coupling with as{M'^) = 0.120 is much bigger than 
the LO coupling with as{Mz) = 0.120, and will do a much better job of fitting the low-Q^ 
structure function data. Hence, the use of the NLO coupling helps alleviate the discrepancy 
between the parton distributions at different orders. Indeed, the NLO coupling is already used 
in some CTEQ LO pdfs and in Monte Carlo generators. 

However, even with this modification the LO fit is still poor compared with NLO, and 
the partons are still depleted in some regions. The problems caused due to the depletion of 
partons have led to a suggestion [11] that relaxing the momentum sum rule for the input parton 
distributions could make LO partons rather more like NLO partons where they are normally 
too small, while allowing the resulting partons still to be bigger than NLO where necessary, i.e 
the small-x gluon and high-x quarks. 
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Figure 8: Comparison between the LO, LO* and NLO up quarks (left) and between LO, LO* 
and NLO gluon (right). 
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Relaxing the momentum sum rule at input^ and using the NLO definition of the strong 
coupling does dramatically improve the quality of the LO global fit (though a K-factor of 1.3 
is necessary for fixed target Drell-Yan data). The = 3066/2235 for the standard LO fit, 
and becomes = 2691/2235 for the modified fit. The data set fit is much the same as in [8] 
and the heavy fiavour prescription used is that in [12], which is a little different from that for 
previous fits at LO. There is a particularly big improvement in the quality of the fit to HERA 
data. When using the NLO definition of the couphng we obtain as{M'^) — 0.121, which is a 
much more sensible result than in the pure LO fit. The momentum carried by input partons 
goes up to 113%. We denote the partons resulting from this fit as the LO* parton distribution 
functions. 

The LO*, LO and NLO partons are compared in Fig. 8. One can see that the LO* quark 
distribution is bigger than at LO at x < 0.1, and only smaller than NLO for a small region 
centred slightly below x = 0.1. Similarly g{x,Q'^) is significantly bigger at LO* than at LO, 
and much bigger than NLO at small x.^ Hence, the LO* gluon distribution should do better 
than LO for gluon-gluon initiated processes (e.g. Higgs production) where ii'-factors are often 
much greater than unity. 

We can make a simple test of the potential of these improved LO* partons by repeating 
the comparisons of the previous section. Consider Drell-Yan production at the LHC, shown in 
Fig. 9. Indeed, we are nearer to the truth with the LO matrix element and LO* pdfs than with 
either LO or NLO pdfs. Moreover, the shape using the LO* pdfs is of similar quality to that 
using the NLO partons with the LO matrix element. So in this case LO* pdfs and NLO pdfs 
are comparably successful. 

The exercise is also repeated for the charm structure function at HERA, as seen in Fig. 10. 
When using the LO coefficient function the LO* pdfs result is indeed nearest to the truth at 
low scales, being generally a slight improvement on the result using LO pdfs, and clearly much 
better than that using NLO pdfs. One would expect there to be a similar result for all processes 
dependent on the gluon distribution at small x , e.g. hadro-production of c and/or h quarks at 
the LHC. 

These simple examples suggest that the LO* pdfs may well be a valuable tool for use with 
Monte Carlo generators at LO, combining much of the advantage of using the NLO distributions 
while avoiding the major pitfalls. However, the examples so far are rather unsophisticated. 
They are sufficient to show that something can go badly wrong in some cases, but in order to 
determine the best set of pdfs to use it is necessary to work a little harder. We need to examine 
a wide variety of contributing parton distributions, both in type of distribution and range of x. 
Also, the above examples are both completely inclusive, they have not taken into account cuts 

^The evolution automatically conserves momentum. 

^Both LO and LO* gluon distributions are smaller than NLO at very large x. This is because the increase 
in the LO quark distributions in this region allows a good fit to Tevatron jet data without a significant high-x 
gluon distribution. The LO gluon distribution at very high x would be more similar to that at NLO in the DIS 
scheme [8]. 
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Figure 9: Comparison of boson production at the LHC using combinations of different orders 
of parton distributions and hard partonic cross-sections, now including LO* pdfs. 
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Figure 10: A comparison of charm production at HERA using LO, NLO, and LO* parton 
distributions with LO coefficient functions and the truth of NLO pdfs with NLO coefficient 
functions. 
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Table 1: A list of investigated processes at the LHC, along with Monte-Carlo codes and the 
hard QCD scales used in the calculations. 

on the data. Nor have they taken account of any of the possible effects of parton showering, 
which is, of course, one of the most important features of Monte Carlo generators. Hence, 
before drawing any firm conclusions we will make a wide variety of comparisons for different 
types of possible process at the LHC, using Monte Carlo generators to produce the details of 
the final state. 

4 Details of Different Production Mechanisms 

We consider a variety of final states for proton-proton coUisions at LHC energies. In each case 
we compare the cross-section with LO matrix elements and full parton showering for the three 
cases of LO, LO* and NLO parton distributions. We also include the results using a Monte Carlo 
generator with NLO matrix element corrections [13] together with NLO parton distributions, 
which we take to be the truth. If a process is not available in MC@NLO, we use the standard 
NLO approximation without the parton shower influence. The different examples span a range 
of values of X and hard scales and probe different combinations of parton distributions. 

Table 1 reports generators ([13], [14], [15], [16], [17], and [18]) used and QCD scales applied 
in the calculations. As the main LO generator we use CompHEP interfaced with FORTRAN 
HERWIG, except pp hh., where HERWIG is more appropriate generator. The main NLO 
generator is MC@NLO except three processes pp Hqq., pp — > Wj — > fivj and pp — > jj, 
which are not available in MC@NLO. The following parameters are used in the calculations: 
mw = 80.4 GeV, mz = 91.176 GeV, Vw = 2.028 GeV, Vz = 2.44 GeV, cos^w- = 0.2311, 
mij — 4.85 GeV, rric — 1.65 GeV, = 105.66 MeV, a — 1/127.9 (fixed), running ag (according 
chosen pdfs). 
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Table 2: The total cross-sections a{pp — > — > jiVn) at the LHC. 
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Table 3: The total cross-sections a{pp — > Wj — > iiv^xj) at the LHC. Cuts: prii) > 20 GeV, 

Wi)l<5. 

4.1 W Production at the LHC. 

We first revisit our previous LHC example of W production, but with full parton showering 
and an investigation of the total cross-section and pt distribution as well as the pseudo-rapidity 
distribution. Hence, we consider the LO production mechanism of quark- antiquark annihilation 
where x — 0.006 at central pseudo-rapidity. The total cross-sections are shown in Table 2, and 
we see that all predictions using the LO matrix elements are lower than the truth, but that 
using the LO* pdfs is easily closest. 

The distributions are seen in the upper of Fig. 11, where the small plots are the ratio of 

the truth to each of our LO matrix element predictions. For the W-boson pseudo-rapidity 
distribution the LO partons give a result which is suppressed at central rapidities due to the 
depletion of quarks seen before. The LO* partons give a good representation of the truth, 
being just a little small at central rapidity and a bit large at the very highest rapidity. The 
NLO partons are slightly worse in normalization, but overall are actually best in shape. They 
are a little low for very high rapidity, where the excess of LO* partons at very high and low 
X overcompensates for enhancement in the NLO matrix elements while the NLO partons do 
not compensate at all. These high rapidities are outside the range of ATLAS and CMS, but 
may be relevant at LHCb. The LO partons give too low a result at central rapidity due to the 
depletion for x ~ 0.01. 

For the pr distributions none of the pdfs gives a good representation of the shape with pt- 
In particular, they all lead to predictions which are too small with increasing large pt- At LO 
the Pt is simulated by a parton shower algorithm (in our case, by the HERWIG algorithm). By 
definition it underestimates the high-p^ region in comparison with the true NLO correction. 
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Figure 11: The comparison between the competing predictions for the differential cross- 
section for VT-boson production at the LHC (left) and for the resulting muon (right). 
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Figure 12: The ranges in Xi and X2 of the contributing parton distributions for VT-boson 
production at the LHC in the different types of calculation. 



As we will see, this type of effect is a recurring feature. It cannot be solved by defining pdfs 
and must be borne in mind whenever using LO Monte Carlo generators. Nevertheless, all pdfs 
produce rather similar shapes using the LO generator, but again the LO* pdfs give the best 
normalization. 

In the lower of Fig. 11 we see the distributions for the final-state muon in the W decay. 
The conclusions are similar. Again for the rapidity distribution the LO* distributions are best 
in normalization, and indeed are an excellent approximation to the truth. The rapidity of the 
muon is usually a bit higher than that of the parent W so the effect right at the edge of the 
rapidity plot in the W case is not visible for the muon, and would be smeared at even higher 
muon rapidity. The NLO partons are low in normalization, but arguably best in shape. Again 
the LO pdfs fail worst at central rapidity. For the pt distribution the NLO excess only sets in 
beyond ~ 50GeV, but this is simply due to the intrinsic pt of the muon from the W decay. 

The range of x and "average" logx = 0.5{logioXi + logioX2) of the contributing partons are 
shown in Fig. 12. For the LO production mechanism x is centred on the line X1X2 = rn^/s 
with some width, whereas at NLO there is a slightly increased contribution from the region 
X1X2 > mvi//y^, where one of the incoming partons emits a relatively hard parton before 
undergoing the annihilation, this hard parton providing some hard pt to balance the W boson. 
A depletion is seen for the LO partons in the region x = 0.001 — 0, 1. 

Since the LHC initial state is dominated by gluons, processes with gluons in the initial state 
are of particular interest. In the W-production process both initial partons are quarks. Initial 
gluons appear in the associated production of W and one jet, and in this process we probe both 
quark and gluon pdfs. We calculated the total cross-section for all our previously considered 



20 



pdf type 


matrix 

plpTTipnf 


(7 \Jlv) ) 


iv-iactor 


NLO 


NLO 


2.40 




LO 


LO 


L85 


1.30 


NLO 


LO 


1.98 


1.26 


LO* 


LO 


2.19 


1.09 



Table 4: The total cross-sections a{pp — > Z/7 — > jiji) at the LHC with cuts {prifJ') > lOGeV, 
\r)iJ,\ < 5.0). 

cases, the LO matrix element with three pdfs and the truth. The NLO approximation was 
calculated by the MCFM program [19]. Reasonable cuts were applied in order to exclude 
kinematic regions where the LO matrix element has singularities. The total cross-sections are 
reported in Table. 3. As we can see, the LO* pdfs give the best normahzation. 

4.2 Z/7 Production at the LHC. 

Next we consider the rather similar process of quark-antiquark annihilation to Z/7 bosons, 
decaying to muons. In order to exclude the dangerous region mi„^ — > 0, where the matrix 
element at LO has a singularity, we apply some experimentally reasonable cuts cuts pr > lOGeV 
and I77I < 5.0. These cuts arc more or less appropriate for most analyses in CMS/ATLAS. The 
process is dominated by the Z peak. The mechanism is rather similar to that for W production, 
but now the quark and antiquark are the same flavour and the x at zero rapidity is slightly 
higher, i.e. Xq = 0.0065. 

The similarity is confirmed by the results. Again all the total cross-sections using the LO 
generators are lower than the truth, as seen in Table 4, but that using the LO* partons is 
easily closest. A similar conclusion holds for the distributions in terms of the final state boson 
or the highest-pr muon shown in the upper and lower plots of Fig. 13 respectively. For the 
boson the LO* partons gives comparable, perhaps marginally better, quality of shapes as the 
NLO partons, but better normalization. The LO partons have the worst suppression at central 
rapidity, and all partons give an underestimate of the high-p^ tail. Similarly, for the muon 
the LO* partons give an excellent result for the rapidity distribution until \r)\ > 4, better in 
shape and normalization that the NLO partons whilst the LO partons struggle at central 77. 
Again, as in W production, the Pt distribution of the muon is better than for the boson, and 
in normalization is best described by the LO* pdfs. 

We illustrate the ranges of Xi and X2 in in Fig. 14. The upper right-hand plot shows that the 
peak in the partons contributing is indeed at X1X2 = m'^z/s, with some contribution from the 
region of higher XiX2. Again, this is particularly for the case of the NLO matrix element, where 
one parton can make a fairly hard emission before the annihilation process. This contribution 
from higher X1X2 is most clearly seen in the lower right-hand plot. The upper left-hand plot 
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Figure 13: The comparison between the competing predictions for the differential cross-section 
for Z/7-boson production at the LHC (upper plots) and for the resulting highest pt muon (lower 
plots) . 
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Figure 14: The distributions of a;i 2 of the contributing parton distributions for Z/7-boson 
production at the LHC in the different types of calculation. 
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Table 5: The total cross-sections for the single top production (with top decay t — > iWrnui b) at 
the LHC. 

shows the contributing x of each parton and is dominated by the line X1X2 ~ Tn%/s, with the 
plateau falling off when the quark distribution for the higher x quark falls away for x — > 1, i.e. 
the contribution is small when xi ^ 0.3 so X2 ~ 0.0001. Perhaps most interesting is the lower 
left-hand plot which shows the ratio of the contribution from particular values of x for the truth 
divided by each of the LO generator calculations. In all cases the ratio diverges at very high 
X due to the possibility of radiation of a hard parton before annihilation at NLO. However, 
over the main region of contribution to the cross-section, i.e. 0.001 < a; < 0.1, the K-factor for 
the NLO and LO* pdfs is flat and close to 1 (though closer for the LO* pdfs), whereas there 
is a distinct bump to values greater than 1 for the LO pdfs. This illustrates very clearly the 
depletion of the quarks in this region for the LO partons, and is very similar in the previous 
case of W production. 

4.3 Single Top Production at the LHC. 

Now we consider a somewhat different process, i.e. single top production with the top decaying 
t ^ /i'^ + i'i^ + b. At the partonic level the dominant interaction process is qb qt (or qb ^ qt), 
where the 6-quark has been emitted from the gluon. Since the 6-quark pdf is determined by 
evolution from the gluon pdfs, this cross-section probes both the gluon distribution and the 
quark distributions for invariant masses of above about 200GcV, i.e. at central rapidity Xq ~ 
0.01 — 0.1. The t-channel nature of this process makes the invariant mass of the flnal state and 
the correspondingly probed x values less precise than the previous case. The total cross-section 
for the various methods of calculation is seen in Table 5 (we excluded Br(t /j,'^ + i'n + b)). In 
this case the result using the LO generators and the LO pdfs is suppressed, but that using the 
LO* pdfs is now larger than the truth. This is due to the large enhancement of the LO* gluon 
distribution. The NLO pdfs give the closest normahzation. 

The distributions in terms of pr and 77 of the final state top and /i originated from the 
top are shown in Fig. 15. For the top distribution the result using the LO generator and the 
LO* and NLO pdfs give a very similar result, being better than the LO pdf result both for 
normalization and for shape due to the suppression of the LO quarks at central rapidities. In 
the case of the fi the distributions calculated with the LO generator look better than for the 
top, since the real NLO correction (irradiation of an extra parton) plays a lesser role for the 
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Figure 15: The comparison between the competing predictions for the differential cross-section 
for single top production at the LHC (upper plots) and for the resulting pt muon (lower plots). 
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Figure 16: The distributions of Xi 2 of the contributing parton distributions for the single top 
production at the LHC in the different types of calculation. 



top decay products. In this process there is a particular NLO enhancement at central rapidity, 
so it gives a total cross-section larger than the truth. In Fig. 16 we display the ranges of x 
sampled. The right-hand plot shows the main difference in this process as compared with the 
previous cases. Here we have much larger values of "average" x. 

The process also illustrates the importance of the NLO corrections in experimental analyses. 
If we accept the picture where qb —>■ qt is the LO approximation for the process^, the main 
correction comes from the gluon splitting to bb, where one 6-quark participates in the top 
production and the second 6 is a parton-spectator. Thus, at LO we have the 6-quark in the 
initial state, but we assume no intrinsic 6-quarks in protons (as we pointed out the b pdfs 
are calculated based on the gluon distribution). It means we must produce the second 6-quark 
somehow in any reliable simulation. Usually Monte-Carlo programs generate the parton via the 
initial parton shower approach. It is clear the matrix element calculation gives a more central 
6-quark with higher pt- So, if the parton is important in an analysis, the LO approximation 
does not work well whatever pdf we apply''. 

^There exists another prescription, where one considers qg qt + b a,s a LO approximation [20]. From our 
point of view, this approach is less motivated both from the mere vertex calculation (both Feynman graphs of 
the process have 3 vertices, whereas qb — s- qt has only two vertices) and availability of the 6-quark pdf. The b 
distribution takes into account some logarithmic corrections which cannot be implemented in the LO matrix 
element. 

""in fact, one can combine the LO approximation and the main real correction due to g ~* bb. See more 
details in [21]. 
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Table 6: The total cross-sections o{jyp — > IT) at the LHC. Strictly speaking this is pp ^ H ^ tt 
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Table 7: The total cross-sections for the process pp — > Hqq at the LHC. 

4.4 Higgs Production at the LHC — Gluon-Gluon Fusion. 

We now consider a process which is determined entirely by the gluon distribution, i.e. pro- 
duction of Higgs of mass 130GeV from gluon-gluon fusion. So the cross-section depends on 
the gluon distribution with the x at central rapidity being xq ~ 0.01. It is well known that 
there is a very large K-factor, approximately 1.7, for Higgs production via this mechanism, so 
it is no surprise that when using the LO generator the cross-section is suppressed by roughly 
this factor using both the LO and NLO pdfs, whose gluon is of a similar size for of a; ~ 0.01. 
However, from the right-hand side of Fig. 8 we see that the LO* gluon distribution is enhanced 
by a factor of ^ 1.25 for the relevant value of x and the extra gluon contribution factor of 
1.25^ compensates a large part of the NLO K-factor. Hence, the result using the LO* pdfs is 
much better than for the LO and NLO pdfs, as seen in Table 6. Since gg ^ H is an s-channel 
process, the "average" x has the same profile as the W and Z/^y production, as we see on the 
right-hand plot in Fig. 18. The distributions in the final state are shown in Fig. 17. The shapes 
are good in all cases, perhaps best for NLO pdfs, but the normalization is poor except for the 
LO* pdfs. 

4.5 Higgs Production at the LHC - Vector Boson Fusion. 

We again consider Higgs production, but via a different mechanism, i.e. a quark from each 
proton emits a vector boson which fuse to produce the final state Higgs boson. In the case 
we use quark pdfs and probe a different x region ^ 0.1. As we can see from Table. 7, in this 
case the NLO K-factor is only a few percent (positive or negative), so the result using the LO 
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Figure 17: The comparison between the competing predictions for the differential cross-section 
for the Higgs production at the LHC (lower plots) and for the resulting r lepton (upper plots). 
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Figure 18: The distributions of X\^2 of the contributing parton distributions for the process 
pp ^ H ^ Tf at the LHC in the different types of calculation. 
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Table 8: The total cross-sections a{pp — > bb) at the LHC. Applied cuts: px > 20GeV, \ri{b)\ < 
5.0, and AR{b,b) > 0.5. 

generator and the NLO pdfs is only slightly above the truth. The result using the LO quarks 
is about 6% too low due to the suppression of the quarks. We note that the LO* quarks are 
similar to the NLO quarks for x = 0.01 so they give a very similar total cross-section, i.e. just 
a little above the truth. The rapidity distributions are very good using both the LO* and NLO 
partons as seen in Fig. 19 but we have, as usual, a small underestimate at central rapidities 
when using the LO pdfs. In this case the NLO corrections do little to the shape of the pt 
distribution. The values of x probed are seen in Fig. 20 and go to quite high values of x. 

4.6 Heavy Quark Production at the LHC. 

We now consider bb and tt production at the LHC. Although the dominant contribution at LO, 
gg/qq ^ bb does not have soft and coUinear singularities (the process is called Flavour Creation 
or FCR), there are two other sources of 6-quarks at hadron colliders at LO: qb —>■ qb, where the 
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Figure 19: The comparison between the competing predictions for the differential cross-section 
for the process pp — > Hqq at the LHC. 



pdf type 


matrix 
element 


a (pb) 


K-factor 


NLO 


NLO 


812.8 




LO 


LO 


561.4 


1.45 


NLO 


LO 


531.0 


1.53 


LO* 


LO 


699.4 


1.16 



Table 9: The total cross-sections a{pp ti) at the LHC. 



30 




Figure 20: The distributions of x\^2 of the contributing parton distributions for the process 
•pj) Hqq at the LHC in the different types of calculation. 



second 6-quark is simulated by initial parton showers (the process is called Flavour Excitation, 
or FEX), and the QCD 2^2 process with light quarks and/or gluons, where the 6-quark 
pair arises in the initial or final parton showers^ (called Gluon Splitting or GSP). These latter 
processes have massless partons and, thus, soft and coUinear singularities. In order to exclude 
the dangerous regions, where the LO approximation does not work, we apply some reasonable 
cuts: prib) > 20GeV, \r]{b)\ < 5.0, AR{b,b) > 0.5. It is interesting that we cannot separate 
the subprocesses at NLO, so only the FOR process exists at NLO (see more details in [22]). 

At the LHC the dominant initial state in this case is gluon-gluon, similar to the inclusive 
Higgs, but at much lower pt (which dominates the total cross-section). So, in the process we 
probe rather low x ~ 10~^ — 10"^ (see Fig. 22). The total cross-sections are shown in Table 8. 
In this case we are at low enough x to be sensitive to the small-x divergence in the NLO 
matrix elements, and the NLO correction is very large. All of the results obtained using the 
LO generator are below the truth, but the reduced NLO gluon means that the NLO pdfs give 
by far the worst result. The best absolute prediction is obtained using the LO* partons. In this 
case the LO gluon distribution is larger than at NLO, so LO pdfs give the second-best result. 

The differential distributions in terms of pt and of a single b quark are shown on the upper 
plots and for the pseudo-rapidity and px of a bb pair on the lower plots in Fig. 21. When using 
the LO generator the LO* pdfs do well for the single b rapidity distribution, but underestimate 
a little at high rapidity. The LO and NLO pdfs are similar in shape, but the normalisation 

^For example, the total cross-section for the improved LO pdfs from Table 8 can be separated into three 
terms: atot — <^fcr. + <^fex + o'gsl, where upcR = 1-6 /^b, apEX — 0.57 /ib, and ctgsp — 0-46 fib - the total 
cross-sections for the FCR, FEX, and GSP processes respectively. 
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Figure 21: Differential cross-sections for b production at the LHC (upper plots) and for a 
pair (lower plots). 
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Figure 22: The distributions of a;i 2 of the contributing parton distributions for the process 
pp bb at the LHC in the different types of calculation. 



is worse. All pdfs obtain roughly the right shape for the ri{bb), except small underestimation 
at very high rapidity. However, for all partons there is a small problem with the shape as a 
function of pt- Obviously, all the ratio curves become higher as pt goes up. As for previous 
processes this happens due to the different behaviour of the additional parton generated at NLO 
compared to those generated by parton showers. In general, we conclude the LO* pdfs give the 
best results in the comparison. Unlike the inclusive Higgs this is not a pure s-channel process, 
and, as we can see in the right-hand plot in Fig. 22, there is no clear peak for "average" x, but 
the main contribution in the process comes from the very low x region, near 10"^'^ — 10~^. In 
this region the LO* approach works better than in the high x region. 

Another interesting heavy quark production process is double top quark production. As a 
short check we calculated the total cross-section for the process. Table 9 reports the numbers. 
At the LHC this process is dominated by the gluon contribution gg tt. For example, 
(^MEiLoi-PDFiLO] = (Jgg^tt + (^qq^tt = 486.9 pb + 74.5 pb. The LO* pdfs appreciably enlarge the 
gluonic cross-section, namely, aME[LO]-PDF[LO*\ = CTgg^tt + CTqq^tt = 622.1 pb + 77.3 pb. Again 
the LO* pdfs give the best prediction. 

4.7 Jet Production at the LHC. 

As our final example we look at high-E'T^ jet production at the LHC. In this case as we span 
the range of we span a range of x, i.e. at LO x = Et/^/s. We also change the dominant 
mechanism as we change E^ - at the highest E^ quark-quark interactions are dominant while 
for a; < 0.1 gluon-gluon interactions take over. In the intermediate range there is also a large 
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Figure 23: Inclusive di-jet cross-section at the LHC Differential cross-sections for the highest- 
Pt jet in the inclusive jet production at the LHC as a function of px (upper plots) and for rj 
(lower plots). 
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Table 10: The total cross-sections cr^pp — > jj) at the LHC with cuts (ptU) > 20GeV, \ri{j)\ < 
5.0, Ai?(j,j) >0.5). 




Figure 24: The distributions of xi^2 of the contributing parton distributions for the inclusive 
2-jets production at the LHC in the different types of calculation. 



contribution from quark-gluon interactions. The jet production as a function of Ex using the 
LO generator and the different pdfs as a ratio to the truth are shown for the lower Et values 
in Fig. 23. For the lowest pt, where we probe x ~ 0.005, the problems due to the relatively 
suppressed NLO gluon are again apparent and the LO and LO* pdfs are better in normalization 
and shape than NLO. Again the LO* pdfs provide the best description, very near to the full 
NLO prediction as seen in table 10. As we go to higher E^ the predictions largely converge, 
but the NLO pdfs used in the LO calculation tend to be a little small. 



5 Discussion 



There are a variety of reasons why the NLO cross-section corrections may be fairly large. As 
discussed earlier in the article, one major reason is when one probes small-a; parton distri- 
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butions and there appears a l/x divergence in the cross-section for the first time at NLO. 
This is particularly the case for gluon dominated processes, and is the main reason for the 
large corrections to the bb production at the LHC, but it also contributes to the lowest Et 
jet cross-section and even probably a little to the Higgs cross-section from gluon-gluon fusion, 
since for the lowish mass used gluons in the region x < 0.001 are probed. There are also large 
corrections due to soft gluon emission, which is important in regions near the edge of phase 
space, i.e. we have large so-called threshold corrections. These have been shown to have effects 
which persist some distance from the strict threshold region, and also contribute to the large 
Higgs production correction in gluon-gluon fusion. They are the dominant reason for the large 
correction in tt production. There can also be a large term appearing in NLO cross-sections 
from the analytic continuation from the space-like region where structure function coefficient 
functions are calculated to the time-like region of s-channel processes [23], and this contributes 
significantly to the large K-factors for W and Z production. 

Hence, even if one extracted the NLO parton distributions perfectly from comparison to 
mainly DIS data, using them along with LO cross-sections for processes at the LHC one would 
expect the predictions to be generally too low, often significantly. This is indeed seen to be 
the case. The particular problem at small x, where the NLO gluon distribution is far too 
small to use with LO cross-sections, was the main motivation for introducing the LO* parton 
distributions, and it is no surprise that they give a far better prediction for the small-x sensitive 
quantities, though the high quality of the quantitative prediction is surprising. However, the 
relaxation of the momentum sum rule has allowed the LO* gluon distribution to be bigger 
than that at NLO and standard LO at most values of x, and the quarks at x < 0.1 driven 
by gluon evolution to follow suit. This actually results in the LO* parton distributions giving 
the best predictions for the quantities where NLO threshold corrections are large, but where 
partons, particularly the gluon distribution, are probed for x < 0.1. In fact, since the parton 
distributions fall away very quickly above x — 0.1 at high scales, the region where x ~ 0.1 
is probed is effectively near threshold. Also, the fact that the increased LO* gluon allows 
the quarks to be bigger than even NLO for x < 0.01 improves the prediction for the vector 
boson production (particularly when an additional jet is present), compensating for the NLO 
correction due to analytic continuation. Overall the LO* parton distributions can qualitatively 
mimic the effect of the NLO cross-section corrections for a surprisingly wide variety of processes. 

However, t-channel dominated processes, illustrated here by single top production and Higgs 
production via vector-boson fusion, generally do not correspond to any of the cases where the 
NLO correction is expected to be large. The term from analytic continuation is not present. 
Moreover, a sharp threshold is not so well defined in a t-channel process, so one does not see large 
threshold corrections in the same way. However, the fact that there is more activity in the final 
state in a f-channel process than a similar s-channel process means that the "average" x probed 
is higher, e.g. nearly an order of magnitude for vector-boson production of Higgs compared 
to gluon-gluon fusion, and there is less likelihood of corrections due to small-a; divergences in 
the cross-section. Hence the absolute NLO cross-section correction is very small in these cases, 
the K-factor when NLO partons are used in both cases being very close to unity, and in fact 
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slightly less. While the LO partons for the relevant x ~ 0.1 are smaller than at NLO, giving 
too small a prediction, the LO* partons are generally enhanced compared to NLO, and give a 
slightly too large prediction. Hence, for our examples of i-channel processes the NLO partons 
are generally the absolute best to use with LO cross- sect ions, but one is only off by at most 10% 
in normalization, with no problem with shape, using either LO or LO* parton distributions. 

6 Conclusions 

We have examined the effects of varying both the order of the matrix elements and the parton 
distributions when calculating cross-section for particle colliders which involve hadrons. The 
intention is to find the best set of parton distributions to use if one only has LO matrix elements, 
as is often the case when using Monte Carlo generators. In order to do this we have calculated 
a variety of processes, both inclusive and more differential, where we can currently find the 
result using both NLO matrix elements and parton distributions. In the case where we look 
at the details of the final state, we have used Monte Carlo generators with either LO or NLO 
matrix elements and LO parton showering. 

We notice that the NLO matrix element corrections are generally positive, and sometimes 
give very large enhancements. In contrast, the parton distributions at NLO are sometimes 
bigger, but also sometimes smaller. Indeed, this must be the case since they satisfy sum rules 
on both parton number and momentum. The fact that higher order partonic matrix elements, 
and sometimes splitting functions contain large terms in ln(l — x) and ln(l/x) means that the 
partons extracted in these regions at LO tend to be enhanced to make up for these missing 
terms. From the sum rules the parton distributions away from these regions are correspondingly 
generally smaller than at higher orders. In practice, there is a pronounced depletion of quarks 
for X ~ 0.1 — 0.001, and the large LO gluon at very small x means it is smaller above x ~ 0.01. 
These are often the regions of most interest at hadron colliders. 

A fixed prescription of either LO or NLO pdfs with LO matrix elements is unsuccessful, 
with each significantly out in some cases. For LO pdfs this is mainly due to the depletion 
just discussed, while for NLO partons the smallness in some regions compared to LO pdfs is a 
major problem if the large NLO matrix element is absent. This is particularly the case at very 
small X, where the large LO gluon compensates for the missing large NLO corrections. The 
best features of each order are required. 

To this end we have suggested an optimal set of partons for Monte Carlos, which is essentially 
LO but with modifications to make results more NLO-hke, and are called LO* pdfs. The NLO 
couphng is used, which is larger at low scales, and helps give a good fit to the data used 
when extracting partons from a global fit. The momentum sum rule is also relaxed for the 
input parton distributions. This allows LO pdfs to be large where it is required for them to 
compensate for missing higher order corrections, but not correspondingly depleted elsewhere. 
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We have compared the LO, NLO and LO* pdfs in LO calculations to the truth, i.e. full 
NLO, for a wide variety of processes which probe different types of pdf, ranges of x and QCD 
scales. In general, the results are very positive. The LO* pdfs nearly always provide the best 
description compared to the truth, especially for the s-channel processes. This is particularly 
the case in terms of the normalization, but the shape is usually at least as good as - and some- 
times much better than - when using NLO pdfs. It is noticeable that the type of enhancement 
of the pdfs provided by the allowed momentum violation seems to give a reasonably universally 
correct modification of the LO pdfs. Hence, while being rather crude, it seems to produce 
the effects obtained by modifying pdfs in a process dependent fashion for use in Monte Carlo 
generators/resummations discussed in e.g. [24] (based on the calculation of the hard Wilson 
coefficient functions for the Qr-divergent part of the cross-sections in the context of Qt resum- 
mation [25]). In particular it produces a relatively small enhancement of quarks and a rather 
larger enhancement of the gluon distribution. 

It should be noted that no modification of the pdfs can hope to reproduce successfully all 
the features of genuine NLO corrections. In particular we noticed the recurrent feature that 
the high-pT distributions are underestimated using the LO generators, and this can only be 
corrected by the inclusion of the emission of a relatively hard additional parton which occurs in 
the NLO matrix element correction. However, if one is limited to a strictly LO calculation, the 
shape as a function of pt seems generally fairly independent of the pdf, and the normalization 
is usually best with the LO* pdfs. 

A prehminary version of the LO* pdfs, based on fitting the same data as in [8], is available 
on request from the authors. A more up-to-date version, based on a fit to all recent data, and 
with uncertainty bands for the pdfs, will be provided in future global fits, i.e. will follow similar 
lines to the MSTW07 NLO pdfs discussed in [26]. 
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